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Abstract. The classical Fibonacci sequence is known to exhibit many fascinating proper- 
ties. In this paper, we explore the Fibonacci sequence and integer sequences generated by 
second order linear recurrence relations with positive integer coefficients from the point of 
view of probability distributions that they induce. We obtain the generalizations of some 
of the known limiting properties of these probability distributions and present certain opti- 
mal properties of the classical Fibonacci sequence in this context. In addition, we also look 
at the self linear convolution of linear recurrence relations with positive integer coefficients. 
Analysis of self linear convolution is focused towards locating the maximum in the resulting 
sequence. This analysis, also highlights the influence that the largest positive real root, of the 
"characteristic equation" of the linear recurrence relations with positive integer coefficients, 
has on the location of the maximum. In particular, when the largest positive real root is 2, 
the location of the maximum is shown to depend on whether the sequence length is odd or 
even. 



1. Introduction 

Starting with any sequence {/[n]}„gN of positive real numbers, we can generate a probability 
distribution on {1,2,..., A^} for every G N as follows: Let X be a random variable taking 
values 1, 2, . . . , such that 



P{X = n) 



f[n] 

N 

Em 

k=l 



In particular, if f[n] are all positive integers, then we can generate such probability distri- 
butions using only integers. In a recent work, Neal [Ij has investigated such probability 
distributions arising out of the Fibonacci sequence and obtained some limiting properties of 
these distributions as — oo. 

While exploring integer sequences from the point of view of generating window functions in 
the context of designing digital filters for signal processing applications [2], we found some 
generalizations of Neal's limiting properties for a class of integer sequences and an optimal 
property of the Fibonacci sequence in this class. In this paper, we present these results. 

2. Probability Distributions induced by Fibonacci Sequence 

The classical Fibonacci sequence [3], [4], is the positive integer sequence f[n] defined by the 
second order recurrence relation 

/[n] = /[n-l]+/[n-2] (2.1) 
1 
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with initial conditions /[I] = 1, f[2] = 1. The ratio of successive terms of the Fibonacci 
sequence converges to the well known golden ratio ip, i.e., 

lim . . = = (2.2 

k^oo J[k\ 2 

For any positive integer N, consider the Fibonacci sequence, /[I], /[2], . . . , f[N], of length A''. 
Using the Fibonacci sequence, we can define a discrete probability distribution as follows. We 
define the probability mass function as 

pM],pJ2],...,pJN] 



Pn(^) = ^^ Vn = l,2,...,iV (2.3) 



IK 

N 

Em 

k=l 

Let X0be a random variable in {1, 2, ... , N} such that 

P{X = n) =p-{n); l<n<N. (2.4) 

In his work, Neal [Ij showed that, E{X), the mean or expectation of X, grows linearly as 
N increases and hence 

lim E{X) = oo (2.5) 

Similarly, let X be a random variable in {1,2, . . . ,N}, where the probabilities, Pjvi'^]' ^o^" 
respond to the Fibonacci sequence, of length A^, in the decreasing order 0. We define the 
probability mass function, 

P(^X = n) = il^j±l^=p^in) ; l<n<N. (2.6) 
k=l 

For such a random variable, 2Li Neal [l| observed the following limiting properties: 

lim E{X) = ^ + l (2.7) 

lim Var{X) = 2(^ + 1 (2.8) 

TV— >-cxD 

where E{X_) denotes the mean or expectation of X and Var(X) denotes the variance of X. 

The aim of this paper, as mentioned earlier, is to show that such limiting properties are 
characteristic for integer sequences defined by linear recurrence relations , of any order, with 
positive integer coefficients. For the sake of simplicity, we shall first discuss the second order 
recurrence relation. 



(.): associated with increasing sequence 

(J : associated with decreasing or the flipped sequence 
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3. Second Order Recurrence Relations 

Let a, b be any two positive integers and consider the integer sequence defined by the second 
order recurrence relation 

f[n]=af[n-l]+bf[n-2] where /[I] = 1; /[2] = 1; (3.1) 

The characteristic equation of this recurrence relation is given by 

- ar - b = 0. (3.2) 

It is well known that the solution to this Eq. ()3.2p is given by 

f[n] = C[Ria,br -Rs{a,br] 

where, 

Ria,b) = 

is the positive real root. The other root of the equation is 



R{a,b) 
and 

Va2 + 46 

We note that, 

R{a, b) > ip for all positive integers o, 6 (3-3) 
R{a,b) = if when a = 6 = l(classical Fibonacci sequence) (3-4) 
We further observe that the ratio of successive terms of this sequence converges to R{a, 6), i.e., 

lim l^ttA = R^a, 6) (3.5) 

fc->-oo J [k\ 

As in section [21 for any positive integer N, consider the above sequence, of length A^, i.e., 

f[l]J[2],...J[N] 
and define probability mass functions, 

as 

PMi^) = ^r^ Vn = l,2,...,iV (3.6) 



Em 



k=l 

Consider a random integer, X, in {1, 2, . . . , N}, such that 

P{X = n)=p-{n) (3.7) 

Similarly let X be a random variable in {1, 2, ... , N}, (with probabilities corresponding to the 
decreasing sequence) such that, 

P^X = n) = ll^:L+l^=p^^n) Vn = l,2,...,iV (3.8) 



Em 



k=l 



In the rest of the paper, for short, we shall denote R{a, 6) by R and Rs{a, 6) by Rs- 
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4. The Main Theorems on the Limiting Properties 
We prove the following limiting properties. Analogous to eq. ()2.5p . we have, 

Theorem 4.1. 

lim E(X) = oo 

N^oo 

Even though E{X) grows as increases, its variance still converges. In fact, we have. 
Theorem 4.2. 

R 



lim Var{X) = 

In the case of A, both E{X) and Var{)C) converge as — >• oo and the limit of Var{)C) 
is the same as that of Var{X). We have, analogous to eq. (|2.7p and eq. (j2.8p . the following limits. 



Theorem 4.3. 

R 



hm E{X) 



N-^oo {R - 1) 

Theorem 4.4. 

R 

lim Var{JC) = — — — ^ 

N^oo {R- 1)2 

We, then, easily observe using these properties that when a = b = 1, the limits of eq. (j4.ip . 
eq. (|4.ip and eq. (|4.ip reduce to the limits, eq. (|2.5p . eq. (|2.7p . eq. (|2.8p . obtained for the Fibonacci 
sequence, by Neal [IJ, described in section [5J Using eq. (j3.3p and eq. (|3.4p . we can prove the 
following optimality properties of the Fibonacci sequence: 

Theorem 4.5. 

max ( lim E{2Q : a, 6 G N ) = lim {E{X) : a = b=l) = cp+l 
Theorem 4.6. 

max I lim Var(X) : a, 6 G N | = max ( lim Var(X) : a,b £ N 



lim {Var{X) : a = b = 1) = 2(p + 1 



We now proceed to prove these theorems. 



5. Some notations used in the following sections 



In order to prove the aforesaid theorems, we make use of certain finite sums, presented in 
Appendix. [A] and Appendix. [Bj We list here, some of the notations used in the subsequent 
sections. 
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Increasing sequence Decreasing Sequence 

N N 



S = Y.f[n] S = Y,f[N + l-n] 

n=l n=l 
N N 



n=l n=l 
N N 

S2 = Y,n^f[n] S, = Y,n'f[N + l-n] 

n=l n=l 

6. Proofs of Theorems related to distributions induced by increasing 

sequences 

We look at the mean and variance of the random variable X, defined as in section [21 
Theorem I4.lt 



Proof. We have, 



lim E(X) = oo 

Af-i>oo 



E{X) = I (6.1) 



Using eq. ijB.SP and eq. ()B.6p . we get 

{NR^+^ -{N + 1)R^+^ + R){R - 1) 



E{X) 



NR- N-1 



R-1 

^E{X) ~ ^-^Zl (6.2) 
E{X) ~ 0{N) 

□ 

Theorem 14.11 now follows from eq. (j6.2p . 
We also observe the following: 

If we now define Y = (N + 1) — X, so that Y is also a random variable taking value 1, 2, 3, 
. . ., N and X has been in a sense centralized, then by eq. (j6.2p . 



R-1 R-1 
Hence 

lim E(Y) = — ^ 

N^oo ^ ^ R-1 

This Y corresponds to X, which we describe in sec. [71 
Theorem 14. 2t 

- R 
lim Var{X) 



N^oc {R-1 



(2 
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Proof. We have, 



Using eq. (fRTO|) and eq- lfRH]) . we get 



E{X') = I (6.3) 



E{X ) (6.4) 



Using eq. (j6.2p and eq. (|6.4p . we get 

Var{X) = E{X^)-E{Xf 



N^{R-lf -2N{R-1)+R + 1 1 



lim Var{X) 



{R-lf \ R-l 

(i? + l) - 1 



Af^oo ^ ' {R- 1)2 

Thus, we get the hmit of the variance as, 

hm Var{X) = (6.5) 

□ 

It is interesting to note from eq. ()6.2p and eq. ()6.5p that, even though the mean increases 
hnearly as the length of the sequence, the variance converges, to a function of R. 

7. Proofs of Theorems related to distributions induced by decreasing 

sequences 

In this section, we look at the probability distribution generated by reversing the sequence 
values. We recall that the probability of the random variable X taking a value n was defined 
as, 

Pi2L = n)= = p^(n) (7.1) 



Theorem 14. 3t 

Proof. We have. 
From eq. fjB.lSp . we get 



Y^f[N + l-k] 

k=l 

1 R 



lim E(X) = 1 + 
A^^oo ^— ^ R-l R-l 

E{X) = ^ (7.2) 

(N + l)S-Si 
S 

~ (AT + 1) - (AT _ -1-) byeq.dSSl) 
hm E{X) = l + -±- = ^ (7.3) 

N^oo K — 1 K — 1 



E{X) 
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□ 

Theorem 14.41 

9, 3R — 1 

lim E{X^) = 1 + 77—^ 

Proof. 

Var{X) = E{X^) - E{Xf (7.4) 
Using eq. ()B.14p we first obtain E{X^). 

§2 



(N + iyS _ 2{N + l)Si S2 

s s s 

{N + if - 2(iV + 1)E{X) + E{X^) 
{R-lf +3R-1 



\2 



lim EiX') = 1 + (7.5) 



From eq. (|7.5p and eq. (|7.3p . we get 

Var{X) = E{X^) - E{Xf 

3R-1 R^ 

1 + 



(i2-i)2 (R-r 



2 



lim Var{X) = ^ (7.6) 

Af->oo [R — \y 

□ 

Thus, we find that, the limit of the variance of the probability distribution with probability 
defined by eq. (j7.ip to be a simple function R, the limit of the ratio of successive elements. 

8. Mean and Variance of discrete probability distributions induced by 

Fibonacci Sequence 

Now, we look at the probability distribution induced by the standard Fibonacci sequence. 
For such a distribution, we obtain the value of mean and variance, using eq. (j6.2p and eq. (|7.5p 
respectively, by using a = 1, 5 = 1, and R = if. 

8.1. Distributions induced by Increasing Sequence. By eq. (|6.2p and eq. (|6.5p . we have 

E{X) ~ 0{N) 

and 

lim Var(X) = , „^ = — ^-ttt 
N^oo ^ ^ (^-1)^ (V^-l)^ 

=> lim Var{X) = 2ip + 1 (8.1) 

Neal [1] looks at the mean and variance of a Fibonacci distribution. It may be pointed out 
that the fact that the variance converges in this case was not observed by Neal. The above 
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derivations show that, in this case, though the mean increases hnearly as the length of the 
sequence, the variance converges to a simple function of the golden ratio ip. 

8.2. Distributions induced by Decreasing Sequence. By eq. (l7.3[) . we have 

lim E(X) = = (/? + ! 

From eq. ()7.6p . we have 

R <f 



lirnVariX) = _ - _ i) 
lim Var{X) = 2ip + I (8.2) 

N^oo 



2 



Comparing eq. (|8.ip and eq. (|8.2p . we find that, in the case of the classical Fibonacci sequence, 
the variance of the distribution converges to the same value, irrespective of whether the prob- 
ability values are increasing or decreasing. 

In the following sections, we make a few observations and remarks related to the classical 
Fibonacci sequence. 

9. Ratio of Variance to Mean 
We observe the curious fact that, from eq. (l4.1j) and eq. M.ip . 



VaAMl ^ _J_ (,,,1) 

N-.00 E{X) R-1 ^ ' 

In the case a = b = 1, that is, in the case of the distribution induced by the Fibonacci sequence, 
since R = ip 

N~.oo E{X) if-l ^ ^ ' 

Thus the ratio of the variance of X to the expectation of X also converges to the Golden 
Ratio, if. 

10. Optimal Probabilistic limit properties of the Fibonacci Sequence 

In this section, we establish certain optimal probabilistic limit properties of the classical 
Fibonacci sequence. 

10.1. Mean. Theorem [OJ 

max ( lim E{X) : a, 6 G N ) = lim (EiX) : a,b = I) = ip + 1 

\N-^oo J N^oo 

Proof. We know that 

R > ip 
=> R-l > ip-1 

11 . . 
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1 



But if - 
Hence, from eq. (|10.ip . we get 



1 



R-1 



^1 + -^ < l + V (10.2) 

< ^ + 1 
R-1 - ^ 

lim {E{X) : a, 6 G N) < lim : a = b = I) (10.3) 

Thus, 

max ( lim E{X) : a, 6 G N ) = lim {E{X) : a, 6 = 1) = 93 + 1 □ 

\N^oo J A"— >oo 

10.2. Variance. Among the family of distributions, induced by second order linear recurrence 
relations with positive integer coefficients, a and b, the Fibonacci sequence gives the maximum 
limit for Var(X), i.e. 
Theorem 14. 6t 

max ( lim Var{X) : a, 6 G N ) = max ( lim VarijQ : a, 6 G N ) (10.4) 

\N^oo J \N^oo J 

= lim iVar{X) : a, 6 = 1) = 2^ + 1 (10.5) 

N—>-oo 

. We prove this as follows: 
Proof. 





R 


> 






1 

R 


< 






1 


> 


1 


- 


R 




=^ 1 - 


1 

R 


> 


1 

1 — 




r 


> 





1X2 



"P. 



From eq. ()10.6p . we get 



<^^ = 2,. + l 



(i?-l)2 - (^-1)2 

R 



(i?-l)2 



< 299 + 1 (10.7) 



lim {Var{X) : a, 6 G N) < lim {Var{X) : a = b = 1) (10.8) 

A"— >oo A— >-cxD 
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From eq. (ll0.8p . we conclude that the limit variance of discrete probability distributions, in- 
duced by decreasing sequence and generated using second order linear recurrence, is maximum 
for the one generated using Fibonacci sequence. □ 

In the next section, we consider the probability distributions generated by higher order 
recurrence relation with positive integer coefficients. 



11. Higher order Recurrence relation with positive integer coefficients 

From Sectionl6]and SectionlTl one can observe that the mean and variance of the probability 
distributions, induced by general second order linear recurrence relation with positive integer 
coefficients, depend only on the largest positive real root, R, of the associated "characteristic 
equation". It turns out that, even in the case of any other higher order recurrence relation 
with positive integer coefficients, it is the largest positive real root that determines the mean 
and variance. We now proceed to prove the same. 
Let 

f[n]=an-if[n-l] + ... + an-kf[n-k] Va,, G N (11.1) 

be a k*"^ order recurrence relation. The characteristic equation of this recurrence relation is 
given by, 

/(x) = x*^ - an-ix''-^ - ... - an-k (11.2) 

Lemma 11.1. Eg. Ul.^) . has atleast one real positive root and if a is the largest positive real 
root, then all the other roots (real and complex) lie within the disc \z\ < a. 

Using Lemma Ill.H it is easy to infer the following: 
With R = a, 

• The same asymptotic relations as in eq. (j7.3p and eq. (j7.6p hold. 

• The same limiting properties as in eq. (j2.7p and eq. (j2.8p hold. 

• The optimality properties of the Fibonacci sequence as in Theorem 14.51 and in Theorem 
14.61 hold good. The reason is that, the largest positive real root, a > (p, the golden 
ratio. 

Thus, Theorem 14.51 can be generalized as 



max lim E{2C) : Oi, z G N = lim (S(X) : ai = 02 = 1, = 0, Vz > 2) = + 1 
and Theorem 14.61 can be further generalized as 
max I lim Var(X) : a,-,i G N | = max ( lim Var(X) : a,-,i G 

\N-^oo ) \N^oo — 

= lim {Var{X) : ai = 02 = l,ai = 0,Vi > 2) = 2^9 + 1 
Now, we proceed to prove Lemma lll.li 
Proof. By Descartes rule of signs, there exists a positive root. Since, 



/(O) < 0, /(oo) > 

Let a be the largest positive root. 
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11.1. a is simple root. Claim: The root is simple. 
Let f{x) be as in eq. (|11.2|) . Since a is a root, we have 

a'' - a„_ia'''"^ + . . . + Un-k = (11.3) 
Suppose, if a is a multiple root, then we must have 

f'ia) = (11.4) 

ka'''^ -{k- l)a„_ia''"2 + . . . + Un-k+i = (11.5) 

Multiplying by a, we get 

ka'^ — {k — l)a„_ia^~^ + . . . + an-k+i<^ = 

ka^ = {k — l)an-ia''^^ + . . . + Un-k+iot 

=^ < k ( an-ia^~^ + . . . + Qn-k+lCt 



=^ < k (^a„„ia'' ^ + . . . + an-k+ia + a„_fc 
=^ a'' < a„_iQ''"^ + . . . + a„_fc+iQ + a„_fc 

This contradicts eq. (jll.3p . Hence a must be a simple root. 

11.2. All roots lie within circle of radius \a\. Now, we prove that all other roots of 
eq. (lll.2p lie within the circle of radius a. Consider eq. (|11.2p . the characteristic equation of 
the corresponding k^^ order recurrence relation with positive integer coefficients. Since a is a 
root, 

/(a) = 

Also, we have 

lim f{x) = oo 

x—^oo 

f{x) > 0, Vx > a (11.6) 

and f{x) > Vx > q (11.7) 
Let /3 be any root (real or complex) of eq. (lll.2p . Clearly, /3 cannot be zero, because, 

/(O) = -an-k < 

Let b=\l3\. 

Claim: b < a 

Reason: Suppose b > a. 

Since /3 is a root, we have 

/(/3) = 0^/3'= = an-i(3^-^ + . . . + a^_k(3^-^ 

= |a„_i/3'=-^ + ... + a„_,,/3"-'=| 

^13" < a„_i|/3|^-i + ... + a„_,|/3r-^ 

^/3''-an-i\/3\''-'-...-an-k\f3r' < 

^ b'' - a^-ib''-^ - . . . - an-kb''-'' < 

^/(6) < (11.8) 

However, since 6 > a, by eq. (|11.7p . f{b) > 0. Thus, we get a contradiction and hence b < a. 
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Claim: b = a 
Since /3 is a root, we have, 

(3^ = an-iP^-^ + ... + a^.fc/S"-'^ 

^1/31" = K-iP"-' + . . . + an-kP'^-'l 

Since we have assumed that \f3\=b=a, 

a'' = |a„_i/3'=-^ + . . . + an-fc/3""''| 

a„_ia'=-^ + . . . + a„_fc = \an-if3'''^ + ■ ■ ■ + an-k\ by eq.([lL3]) (11.9) 
an-i\l3\''-^ + . . . + an-k = |a„_i/3^-i + ...+a„_fc| (11.10) 

=^ a„ — . . . , ttn-k are ah complex numbers in the same direction. 

=^ an-k,CLn-k+i(3 are real positive 
=^> /3 is real positive 
=^ (3 = a. 

Thus, based on the above claims, we find that 

• the only root on the circle \z\ = a is z = q and there is no root outside this circle 

• all roots not equal to a are inside the circle \z\ = a. 

11.3. a is strictly greater than the golden ratio (p. We prove that the largest positive 
real root of the characteristic equation, corresponding to higher order linear recurrence relation 
with integer coefficients, is greater than the golden ratio (p. 

= an-ia''-^ + . . . + an-ka''''' 

Since all terms on the RHS are greater than 



9 

a > a + 1 



As a„_i, a„_2 are positive integers, we see that, 

> (11-11) 

From the fact that, the golden ratio, is the largest positive root of — x — 1 = and that 
x^ — X — 1 > for X > it follows that a is always greater than ip. □ 

12. Convergence in Distribution 

It must be noted that the distribution {p^{n)},l < n < N, converges to a geometric 
distribution on N = {1,2, 3,...}. More precisely, we have the following: 
Let y be a random integer in N, such that 

p(n) = P{Y = n)=p'' (^y^) (12-1) 

where p = We have, 



E/'" = T^ (12.2) 



n=l 
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and hence 



n-l d ( _P_\ 1 



Thus, 



Consequently, 



oo 
n=l 



(1 - Pf 



oo ^ 

E{Y) = Y,np-i 

n=l ^ 



1-p 



P 



(1 - pf 
1 



1-p 



1-p 
1 



1 



Similarly we can see that 
Thus, we see that. 



^E{Y) 
Var{Y^) = 



R 



R-l 



R 



{R-l) 



(12.3) 
(12.4) 



hm E{X) = E{Y) 
lim VariJQ = Var(Y) 

N^oo 

We now show that the random variable X converges to the random variable Y in distribution. 
We have for every n G N, 



\ R^+^ - R 
1-R-^ 



1-p 

p" 



N 



N 



Thus, 



1-p 



1-p^ 



{R-l) 

(R-l) 
1 

- - 1 
P 

P 

1-p 



Since lim p =Oasp= — <1, we have 

N-^oo R 



lim p^{n) = p" 

N^oo — 



n f^-p 



(12.5) 
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From eq. (ll2.ip and eq. (ll2.5|) we see that 

lim 

N-¥oo — 

and thus X converges to Y in distribution. 



lim p^{n) = p{n) Vn G N 

N-¥oo — 



(12.6) 



13. Convergence of Moments 

In the previous subsection, we observed that the distributions generated by second order 
Unear recurrence relation converge to a geometric distribution, i.e., 

• X converges to Y in distribution and 

• E{X_) converges to E{Y) 

• E(X'^) converges to i?(y^) (and hence Var(X) converges to Var{Y)). 

Now, we prove that all the moments of X_ converge to the corresponding moment of y, as 
N ^ oo. 



Proof. For any positive integer k, we have, 



N 



n=l 



N 



n=l 



We have. 



Also, 



N 



n=l 



N+l-n 



R 



N+l 



R 



R-l 



N 



N 



N+l-n 



n=l 



n=l 



Substituting eq. ()13.2p and eq. p3.3p in eq. ()13.ip , we get 



Eix! 



n=l 



- R 



\ 



(R-l) 



J 



f JL \ 

n=l 

R^+^l-R-^) 



\ 

N 
n=l 

1-R-^ 



(R-l) 



J 



1-p 



(13.1) 



(13.2) 



(13.3) 
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Hence, as N ^ oo, 

N /I _ \ 

Eix!") ^nVf— because R-^ ^ 

= E{Y'') 

Thus, we have for every non-negative integer k, 

showing that all the moments of X converge to the corresponding moments of Y. □ 

13.1. Limits of Moments in the case of classical Fibonacci sequence. The moments 
of the probability distribution induced by classical Fibonacci sequence converge to a linear 
function of the Golden ratio, 

Proof. We have, 

oo 

Ex"- = for Ixl < 1. 

1-x 

n=l 

which on differentiation gives 

oo ^ 



n=l 



From eq. ()13.4p . we get 



Then, 



n=l 



oo 



n=l 

Similarly, differentiating eq. (|13.5p . we get 

oo 



n=l 

[l-x) + 2x 1 + X 



"■a; = 7 -7 



and hence, 



and thus, 



(l-x)3 (l-x)3 
(1 + x)x 



n=l 



3 



n=l 



(13.7) 



E"=^" C-^) - ^ (13 ^ 



We shall now show, by induction, that in general, for any positive integer k, 

n=l ^ / \ I 
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where P^_j^{x) is a polynomial in x of degree k — 1. We have from eq. ()13.6p and eq. (ll3.7p . that 
the result eq. ()13.9p is true for k = 1 and k = 2. Thus, it is enough to show that if eq. ()13.9p is 
true for some k = r, then it is also true for k = r + 1. Suppose eq. ()13.9p is true for k = r, then 



n=l 

Hence, we get, 



l-x\ P_,ix) 



:i-x) 



On differentiating eq. (jl3.1Up . we get, 



^ (1 - x)2'"+2 

(1 - x){P^_, (x) + xP;_^ (x)} + xP^_, (x) 



r 



(1 - x)'-+2 



(13.11) 



(1 - x)'-+2 

where 

PXx) = (l-x){P_,(x) + xP'_^(x)} + xP,_,(x)r (13.12) 
is a polynomial of degree r. From eq. (jl3.1ip . we get, 

xPix) 



Z]'" - (1 _ xy+2 

n=l 



n'^+^x" 



and hence, 

'l-x\ PXx) 



1 - X 



ir+l 



n=l 

and thus proving that if eq. (jl3.9p is true for k = r, then it is also true for k = r + 1. Thus by 
induction, we get eq. (|13.9p is true for all k. From eq. (|13.9p . we get 



n=l \ P / 



and hence 

lim E{X^) = EiY") = (13.13) 



N^oo ' ^ ' (1 - p)'' 



eq. (jl3.13p is a rational function of p and hence a rational function of R. In particular, for the 
classical Fibonacci sequence, we have 
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Since — = 99 — 1, we get 

= Qi^p), a polynomial in cp (13.14) 

But, since ip'^ = ip + 1, all powers if^ for j > 2 can be expressed as linear polynomials in ip and 
hence we get any polynomial in ip can be rewritten as a linear polynomial in (p. Thus, 

lim E{X'') = E{Y) (13.15) 

= Ak^ + Bk (13.16) 

□ 

We now provide a detailed analysis about the location of the maximum in the sequence 
resulting from self linear convolution of linear recurrence relations with positive integer se- 
quences. 

14. Linear Convolution of linear recurrence relations with themselves 

In the previous sections, we proved certain optimal properties of the classical Fibonacci 
sequence among the general fc-th order recurrence relations with positive integer coefficients, 
with respect to the discrete probability distributions that they induced. In this section we bring 
out another fascinating property of the classical Fibonacci sequence. Let be a finite length 
integer sequence, of length A^, generated by a k-th order recurrence relation with positive 
integer coefficients. The location of maxima in the sequence, y„ resulting from the linear 
convolution of with itself, can be only one of the following: 

• 2A^-1, when the largest positive real root of the "characteristic equation" is less than 
2 

• 2A^-1, when the largest positive real root of the "characteristic equation" is equal to 2 
and the length even. 

• 2N when the largest positive real root is equal to 2 and the sequence length is odd. 

• 2N when the largest positive real root is greater than 2. 
We prove the above in the following sections. 

15. Some Notations and Definitions 

As mentioned earlier, while investigating the use of integer sequences for window functions 
in digital filtering [2|, we observed that most of the integer integer sequences generated by k- 
th order recurrence relations, with positive integer coefficients, are non-decreasing. However, 
conventional classical window functions are symmetric, bounded and smooth. Though, integer 
sequences, mentioned as above, are bounded for finite length, they are not symmetric. Hence 
to be used as window functions, integer sequences are symmetrized. However, this does not 
guarantee a smooth profile for most of the integer sequences. In order to get a smooth profile, 
we do one of the following: 

• Perform linear convolution after symmetrizing the integer sequence about, say, N/2. 

• Perform linear convolution of the non-decreasing integer sequence and then symmet- 
rically extend the resulting sequence about its maximum. 
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While in the first case, for any finite length, A^, one can symmetrize at half the length, 
the second case requires the location of maximum in the sequence resulting from the linear 
convolution. In the sections to follow, we obtain the location of maxima in the sequence 
resulting from the linear convolution of with itself. In particular, we look at the family of 
sequences generated using fc-th order linear recurrence relation. First, we look at the second 
order case and then generalize the result to higher order recurrence relations. 



15.1. Notations. Let 



N 

X = {Xjn=l (15-1) 



, 2JV 



y = {yJn=2 (15.2) 

x„ = aa;„_i + 6x„_2, Va, 6 G N (15.3) 
The ^^characteristic" equation of the above is given by 

- ax - 6 = (15.4) 
The largest positive real root of the above equation is as below: 



i? = a+^?±5 (15.5) 
15.2. Definitions. The linear convolution of two sequences is given by 

n-l 

yn=E^'=^n-^ 2<n<iV (15.6) 

k=l 

N 

k=j 



For 2 < n < N, we have, 

y^-Vn-i = - Ex,x„_^_, (15.^ 



n— 1 n—2 



k=l k=l 
n-2 

= -2;„_i_J+a:„_iX, (15.9) 
k=i 

> 

> i<n<iv (15.10) 

From eq. ()15.10p , we observe that the convolution sum is an increasing sequence for 1 < n < 

N. 

Now we look at the eq J15.7l for j > 1. 

N N 

yN+j ~ Vn+j-I ~ y~^^fc^Ar+j-fc ~ Xf,Xj^^j_-^_f, (15.11) 

k=j k=j-l 

At this juncture, we revisit the general second order recurrence relation with positive integer 
coefficients, given in eq. (|15.3|) . We evaluate eq. p5.3p for different possible values of a and b 
and hence compute the difference as in eq. (ll5.1ip . 
(1) a = l 

(a) 6 = 1 
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(b) 6 = 2 

(c) 6 > 2 
(2) a > 1 

16. g = 1,6 = 1 : 

We look at case (fTa| that corresponds to a = 6 = 1. Eq. (jl5.3p . in this case, generates the 
elements of classical Fibonacci sequence. Eq. (jl5.3p can be rewritten as 

= + (16.1) 

R, the largest positive real root in the case of classical Fibonacci sequence corresponds to the 

golden ratio, (f = ' With the above conditions, eq JlS.llI gives 

N N 

Vn+j ~ Vn+j-I ~ ^^^fc^JV+j-fc ~ ^k-^N+j-l-k (16.2) 

k=j k=j-l 
N N 

'^^•^ki'^N+j-l-k ~^ •^N+j-2-k) ~ ^ ^ ■^k-^N+j-l-k 
k=j k=j-l 

N 
k=j 



•^k'^N+j-2-k -^j-l-^N 



Suppose, J < — 1, 



Rearranging the terms, we get 

We find that the term on the right hand side, in the above equation is greater than 0. Hence 

Vn+j - Vn+j-i > (16-3) 

From eq. (jl6.3p . we find the sequence to be increasing for 1 < j < — 1. The only pair of 
elements of y,^ that needs to be compared are y^^ and y^M-i- From eq. (|16.2p . we get 
For N — 1 < j < N , the above corresponds to 

y2N-y2N-i = xl-2x^x^_, (16.4) 
= x^{x^ — 2xj^_-^) 

= x^(x^_2 —x^_-^) since x^ =Xj^_-^ ~^x^_^ 
We know that x„ < x^_^_-^ , and hence, from the above equation, we have 

^ - y2N-i < (16.5) 

From eq. (|16.5p . we find that y^^f < y2jv-i- Thus, in the sequence resulting from the linear 
convolution of finite length classical Fibonacci sequence with itself, the maximum is located 
at 2 — 1 , where N is the length of the sequence x^ . 
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17. g = 1,6 = 2 : 

Now, we explore what happens when a = 1,6 = 2 in eq. (|15.3|) . as mentioned in case ijlbp . 
The corresponding second order recurrence relation is as below: 

Xr,= X„_^+2Xn-2 (17.1) 

In this case, the ''^characteristic" equation is j;^ — x — 2 = and the largest positive root, R, 
is 2. To obtain the location of the maximum in the linear convolution of eq ll7.1l with itself, 
we adopt the same procedure as classical Fibonacci sequence. 

N N 

k=j k=j~l 
N N 

k=j k=j-l 
N 

~ ^^2a;^x^^^_2_j. — x-_-^^Xj^ (17.2) 

k=j 

(17.3) 

Suppose, j < N — 1, 

— ~^ -^ivl^j-l ~ ^X^.g) — X^_-^X^ 

= — 2x._^Xj^ 

= 2xj^_-^{x._^ + 2x._^) — 2x._^{x^_-^ + 2xj^_2) 

Rearranging the terms, we get, 

= 2^]v-i(-^j-2 ~ 2a;^_3) + 4x^._3 —x^_^) 
Vn+j - Vn+.-i > (17.4) 
From eq. (jl7.4p . we find the sequence to be an increasing one for 1 < j < N — 1. The only 
elements of y,^ that need to be compared are y^^f and yaiv-i- 
For — 1 < j < N, the above corresponds to 

y2N~y2N-i ~ x^ — 2x^Xj^_-^ (17.5) 
= x^{x^ -2xj^_^) (17.6) 

• When N is even, < 2x^_-^ and hence yj^y-i maximum. 

• When N is odd, > 2x^_-^. Therefore, y^j^ is maximum. 
This can be easily proved by mathematical induction. 

Proof. Let k = 2. With standard initial conditions 1, 1 one can observe that /2 < ^/s and 
/2 < 2/i. 

Let = 4. The sequence elements with initial conditions 1, 1 turn out to be 1, 1, 3, 5. We 

observe that < ^f^ and /4 < 2/3. 

For any k, we prove the following by induction. 

<min{2/,,_,,^/,,^J (17.7) 
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Let eq. (ll7.7p be true for k-1. So, we have 

<mm{2/,,_3,i/,,_J (17.8) 
We now prove that eq. ()17.7p is true for k i.e., 



= A._i+2/,,_, (17.9) 



We know that from eq. (|17.8p . /„.^ „ < -/a^^i 



Similarly, 



Hence/2, < + 

^A. < 2/,,_, (17.10) 



(17.11) 



/ 2fc+l 




f2k + 2/2,_i 


< 2/,,_, 






Hence/2, 


< 


/2fe+l ~ Afc 


^2/2, 


< 


/ 2fe+l 










< 


1/ 

2 ■' 2fc+l 



(17.12) 

Thus we find that, f^^, < min{2/2,_-^, ^f2k+i}- From this one can see that has its maximum 
either at 2A^-1 or 2N depending on whether the sequence length is odd or even. □ 



18. g = 1,6 > 2 : 

With a = 1,6 > 2 in eq. (|15.3p . we look at case (fTc|) . The corresponding second order 
recurrence relation is as below: 

= x„_-, +bXn-2 (18.1) 

In this case, the ^^characteristic" equation is — x — b = and the largest positive root, 
R, greater than 2. We now proceed to obtain the location of the maximum in the linear 
convolution of eq. (jl8.ip with itself. 

N N 

k=j k=j-l 

N N 

k=j k=j-l 

N 

= X]^^fc^iv+i-2-fc -Xj_,x^ (18.2) 

k=j 
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Suppose, J < iV — 1, 

Un+j ~ UN+j-l — ^-^iV-l-^j-l + bx^X._,^ — X._-^Xj^ 

> bx^_,x^_^ + bx^{x._, - bx^_^) - x^^.x^ 
— ^^iv-l-^j-l ~ bx-_^Xj^ 

= tejv_i(2;,_2 + bx^_^) - 6x^_3(x^_, + hx^_^) 



Rearranging the terms, we get, 

= bx^_^{x - X, 3) + 6^x, 3(xj^_i - x^_^) (18.3) 



J- 

> (18.4) 



From eq. (|18.4| ). we find the sequence to be an increasing one for 1 < j < N — 1. The only 
elements of y„ that need to be compared are y^^^ and y^^^i- 
For — 1 < j < N, the above corresponds to 

~ y2N-l = ^If ~ '^^N^N-l (18.5) 

= x^{x^_^+bx^_^-2x^_^) (18.7) 

Let us consider the sequence generated by eq. (|18.ip . with 6 = 3. The elements of the sequence, 
with initial values being 1, 1, are 

1, 1, 4, 7, 19, 40, 97, ... 

From the above, we find that bx^_^ is always greater than x,^_-^ . Thus eq J18.7l is always positive. 
Hence, we find that y^j^ > y2jv-i- Thus, the maximum is located at 2N and not at 2N — 1. A 
formal proof by induction can be obtained along the same lines of mathematical induction, as 
found in section. [T71 



g > 1 : 

We now present the case when o > 1 and 6 G N, in eq. (|15.3p . 

N N 

Vn+j ~ Vn+J-I ~ ^ ' ' j -^k-'^N+i-k ~ ^ > ^k-^N+j-l-k 

k=j k=j-l 

N N 

k=j k=j-l 
N N 

= '^{a-l)x,^x^^._^_^+^bx^x^^^_^_^ -x^_,x^ (18.. 

k=j k=j 



-■^ N+j-\~k 



(18.9) 

N 

When a > 1, the term x._^Xj^ gets accounted for by ^^(o — l)Xj.Xj^^j._-^_j,, with k = N. The 

k=j 

remaining terms are all positive. Hence y„ increases till 2N — 1. Now, we compare yjjv-i 
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y^pf- For N — 1 < j < N, the above corresponds to 



- y2N-i = ^l- '^x^x^_^ (18.10) 

= x^{x^-2x^^,) (18.11) 

= x^{ax^^^+bx^^^-2x^^J (18.12) 

= x^{{a-2)x^_^+bx,^_2) (18.13) 



With a > 1 and 6 G N, eq. (ll8.13p is always positive. Hence y(n) has its maximum at t/jjv- 

19. Self linear convolution of higher order linear recurrence relations with 

integer coefficients 

In the previous section, we found that the sequence resulting from the self linear convolution 
of second order linear recurrence (with positive integer coefficients) has its maximum located 
at either 2N-1 or 2 A'" is determined by the 

• coefficients 

• and hence the largest positive real root. 

It was proved that, when the largest positive real root was 2, the location of maximum is either 
2A'^-1 or 2N, based on whether N, the length of x^^, is even or odd. It is also worth noting that 
this is true even in the case of self convolution of, x^ = ax^_-^, first order linear recurrence 
relation (with positive integer coefficients). In this case, if a = 2, the largest positive real root 
is 2. has its maximum at either 2A'^-1 or 2N depending on whether the length is even or 
odd. However, when a = 1, it can be easily proved that, the maximum of y^ is located at N. 
Now, we look at the higher order recurrence relations with positive integer coefficients. Let 

x„=a,^^iX^^, + ... + a^x^_^ Va, GN (19.1) 

be a higher order recurrence relation. On similar lines, as self linear convolution of the sec- 
ond order linear recurrence relation, one can establish that even in the case of higher order 
recurrence relations, there exist three cases, namely, 

• all the a. 's are 1. y^ is maximum at n = 2N — 1 

• a. 's are 1 and = 2. y^ is maximum at n = 2A^ — 1 or n = 2N , depending on whether 
the length is even or odd 

• o-'s are greater than 1. y„ is maximum at n = 2N 

20. CONCLUSION 

In this work, on the one hand, we have explored and proved certain optimal probabilistic 
limit properties of Fibonacci sequence. We have proved that, among the distributions induced 
by A;*^ order linear recurrence equations with positive integer coefficients, the limits of the 
mean and the variance are maximum for those generated by second order linear recurrence 
relation, corresponding to the Fibonacci sequence. We also found that the ratio of the variance 
to the mean, of the distribution induced by the Fibonacci sequence turns out to be golden 
ratio, ^p. Moreover, for all such distributions, we have, 

• X_ converges to Y in distribution 

• all the moments of X_ converge to the corresponding moments of Y 

• in the classical Fibonacci case, all moments of X converges to linear function of the 
golden ratio, (p. 
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Besides proving certain optimal probabilistic limit properties, we also analyzed the second 
order linear recurrence relation with positive integer coefficients, from the point of view of self 
linear convolution. We established that the sequence resulting from the self linear convolution 
of any higher order linear recurrence relation with integer coefficients, can have its maximum 
either at 2N or 2N-1. A more interesting observation in this regard is that when the largest 
positive real root of such recurrence relations is 2, the location of maximum depends on whether 
the sequence length of x„ is odd or even. 

Appendix A. Some Preliminaries 
We list some of the identities used in this paper in the following sections. 



N 



G{x) = ^x" = x + x^ + ... + x^= (A.l) 

n=l ^ ^ 

Let G'{x) be the derivative of G{x). We have, 

jv-i _ (x-l){(JV + l)x^-l}-(x^+i-x) 



G'{x) ^l + 2x + ... + Nx 



Hence, we get. 



TV 



{x-lf) 

{N + l)x^+^ ^x^{N + l)x^ + 1 - x^+i + X 

{x-iy 

(X- 1)2 



a;G'(x) = ^nx" = x + 2x^ + . . . + Nx'^ 

-{N + l)x^+i + 



n=l 



(A.2) 



From this we get. 



(X- 1)2 

{xG'{x))' = l + 2^x + ... + N^x^~^ (A.3) 



which by eq. (|A.2p is 

{x - lf{N(N + 2)1^+1 -{N + ifx^ + 1} - 2{x - l){Afx^+2 - {N + l)x^+^ + x} 



{xG'{x))' 



{x - 1)4 

(x - l){Ar(Ar + 2)x^+i - {N + l)^x'^ + 1} - 2{Afx^+2 - {N + l)x^+'^ + x} 

{x - 1)3 

N2^^n+2 _ 23.JV+1 _^ ^JV} _^ N{-2x^+'^ +2x^} + {x^+^ +x^ + X + 1} 



{x - 1)3 



and hence, 

N 

x{xG'{x)y = 



n=l 

jV2{^iV+3 _ 2-^iv+2 + 3;iv+i} + Af{-2x^+2 + 2x^+1} + + x^+i + x^ + x} 

(x - 1)3 

Appendix B. Useful Finite Sums 



{A.4) 



In this section, we obtain the expressions for some of the finite sums relevant to the work 
presented in this paper. 
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N 



B.l. f[n]. Let S be the sum of the sequence f[n] defined as 



n=l 

N N 



n = l 71: 

From eq. ()A.ip . we get 



S = Y f[n] = y , ^ :(fi"-i?g) (B.l) 



= i? + i?' + i?^ + ... + i?^ = ^i?"= (B.2) 



R-\ 



^ <+i - Rs 



G{Rs) = Rs + Rl + Rl + ... + R^ = Y.Rs= \ _i (B.3) 

n=l 

Substituting ea. dR^ and ea. dRSll ea. dBl^ . we get, 

1 ^ 

1 ( r^+^-R R^+^-Rs 
V^^TM \ R-l Rs-l 

We have, 



-R = : 



R^ > b 



B.2. ^^n/[n]. Let Si be defined as 

n=l 

N 

Si = ^nf[n 

n=l 

N 



^ n=l 



(BA) 



Hence, taking i?^^^ common from eq. ()B.4p . we see that, 

Va2 + 46 \R-lJ 

N 



By eq. (IA.2j) . we get 

f NR^^'-{N + 1)R^^^ + R 

l^nR ^^^-^p (B.7) 

n=l 



26 ARULALAN RAJAN, R. VITTAL RAO, ASHOK RAO, AND H. S. JAMADAGNI 

We can get a similar expression involving the other root Rs- Using the fact that the absolute 
value of the ratio of Rs and R is less than 1, we obtain 

1 f NR^+^-{N + l)R^+' + R \ 

^ 7m I ^TTTp ) (^-^^ 

TV 



B.3. ^^n^/[n]. Let ^2 be defined as 

n=l 

N 

S2=^n^f[n] (B.9) 



n=l 

Using eq. ()A.4p we get 

N^rN+s _ 2RN+^ + + N{-2R^+^ + 2i?^+i} + {R^+^ + R^+^ + R"^ + R} 



So 



{R-lf 



(B.IO) 



N 



B.4. ^ /[iV + 1 - n]. Let S be defined as 

n=l 

TV 

S = J]]/[iV + l-n] (B.ll) 

n=l 

N N 

Clearly, S = ^f[n] = ^/[iV + l-n] = 5 (B.12) 

n=l n=l 
Thus the sum S_ = S , where S is as in eq. (IB.6P . 



TV TV 

n] 

n=l n=l 



B.5. nf[N + 1 - n]. Let denote ^ n/[iV + 1 

S, = Y,nf[N + l 

n=l 

TV 
n=l 

TV 

- J];(iV + 1 - n)f[N +l-n] + {N+ 1)^ /[iV + 1 



TV 

n] 

n=l 

TV 
n=l 

TV TV 

n=l n=l 

TV 



-J^fe/[A;] + (iV + l)5 

k=l 

>S^ = {N + l)S-Si (B.13) 
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N 

B.6. '^n'^flN + l-n]. Let S2 be defined as 

n=l 

N 

71=1 

N 

= Y^{N+l-nff[n] 

n=l 

N N N 

= {N+ 1)2^ f[n] - 2{N + 1)^ nf[n] + ^ n'f[n] 

n=l n=l n=l 

^ 3^ = {N + 1)^S -2{N + l)Si + S2 (B.14) 
where 5i and S2 are as in eq. flB.Th and in eq. flB.Qh respectively. 
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